Quantum capacitance: a microscopic derivation 
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We start from microscopic approach to many body physics and show the analytical steps and 
approximations required to arrive at the concept of quantum capacitance. These approximations 
are valid only in the semi-classical limit and the quantum capacitance in that case is determined 
by Lindhard function. The effective capacitance is the geometrical capacitance and the quantum 
capacitance in series, and this too is established starting from a microscopic theory. 
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I. INTRODUCTION 

Several new concepts and ideas have developed in last 
few decades on nano-electronics and they are often ques- 
tioned [l| . AC response of quantum dots in the coherent 
regime has been measured in recent experiments • A 
good understanding and control over such phenomenon 
can lead to many novel devices, specially in metrology 
0, 0] ■ The experimental results are analyzed in a series 
of works, using effective variables like quantum capaci- 
tance |7H12| . Capacitance of mesoscopic systems are very 
different from geometric capacitance. In mesoscopic sys- 
tems one can differentiate between electrostatic capac- 
itance and electro-chemical capacitance. Although, in 
principle one can also define a magnetic field induced ca- 
pacitance, in practice one defines a field dependent elec- 
trostatic or electrochemical capacitance [l3j]. Ref [l4| 
gives a detailed analysis of electrochemical capacitance 
which gives corrections to the geometrical capacitance 
due to field penetration into the conductor which occurs 
over a finite length scale comparable to the dimensions of 
the sample and ignored in large systems. Electrochem- 
ical capacitance is a property of open systems (systems 
connected to leads and electron reservoirs). In such open 
systems, electron-electron interaction cannot be treated 
exactly and characteristic potentials were introduced to 
account for Coulomb interaction approximately. The cor- 
rection term appear as another capacitance in series with 
the geometric capacitance. Both, open and closed sys- 
tems can have an electro-static voltage induced capac- 
itance. A closed system of a finite wire (referred as a 
stub) connected to a closed ring was analyzed in Ref. 
[l3|. The system was reduced to a two level system, 
wherein there is a hybridization of a single level coming 
from the stub and another single level coming from the 
ring. Coulomb interaction was again treated approxi- 
mately with the help of characteristic potentials and sin- 
gle particle level approximations. Quantum corrections 
was again shown to appear as a capacitance that appear 
in series with the geometric capacitance. The quantum 
capacitance is given by the Lindhard function [13| . Sub- 
sequently, several authors have tried to interpret experi- 
mental data and numerical calculations in terms of quan- 



tum capacitance (l5l . Il6j . A microscopic analysis stating 
under what circumstances and assumptions one can use 
such a parameter is not done so far. 



Capacitance of a system is self consistently determined 
by Coulomb interactions and this is no exception for 
quantum capacitance as well. However, quantum me- 
chanically electrons can also interact via Fermi statistics 
and so even when Coulomb interaction is ignored, a sys- 
tem can have a quantum capacitance. While geometric or 
classical capacitance is determined by the volume, shape 
and dielectric constant of the system, charge in quantum 
mechanics can reside in orbitals that do not have a space- 
time description. The existence of an effective variable 
of quantum capacitance, can simplify the complexity of 
many body physics. 



Unlike that in open systems, electron-electron interac- 
tion and statistics can be treated exactly in closed sys- 
tems. In this work we deal with closed systems so that 
an analytical proof can be given and exact numerical di- 
agonalization is possible for verification. We would like 
to analyze the assumptions and concepts required to ar- 
rive at a statistical mechanical variable of electrostatic 
quantum capacitance. When a system is weakly cou- 
pled to a reservoir, making it an open system, one can 
describe the system in terms of the eigenenergies of the 
system and Fermi-Dirac distribution function [l7| . So 
our results are also valid for weakly coupled open sys- 
tems. Electro-chemical potential also works by affecting 
the electrostatic potential inside the system So if 

electrostatic quantum capacitance cannot be proved then 
electro-chemical quantum capacitance may also not hold. 



Our analysis is independent of model and valid for any 
arbitrary geometry in any dimension. However, we will 
use some models and systems for numerical verification, 
that are described in section II. Analytical derivation of 
quantum capacitance is given in section III. Conclusions 
are given in section IV. 
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II. MODELS FOR NUMERICAL 
VERIFICATION AND ILLUSTRATION 

We have given in Fig. 1, schematic diagrams of a ring 
(Fig. la), a stub connected to a ring (Fig. lb), and a 
2D square geometry (Fig. lc). Although our analysis is 
not restricted to these geometries, we will use them as 
reference and examples. Figs, (la), (lb) and (lc) repre- 
sents continuum cases whereas Figs. (Id), (le) and (If) 
represent discrete versions of the same systems as that 
in (la), (lb) and (lc), respectively. Discrete models are 
useful for numerical analysis. The vector potential due to 
a magnetic field can non-trivially change the electronic 
states of a system (due to quantum interference) while 
having very little effect on the bound positive charges 
that can be assumed to be uniform [13j . It is very easy 
to see polarization due to vector potential in rings as the 
magnetic field can remain confined to the center of the 
ring while the electrons in the ring feel only the vector po- 
tential. However, it also occurs in small two-dimensional 
or three-dimensional quantum systems where weak mag- 
netic fields have negligible effect (Lorentz's force being 
weighted down by the velocity of light), but vector po- 
tential will drastically change the states of the system 
due to interference effects. For numerical verification, we 
use the generalized Hubbard Hamiltonian describing a 
discrete system consisting of sites. 
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where ti is the site energy of the ith site, t is the hopping 
parameter (in presence of magnetic field it becomes com- 
plex, i.e., t — > texp[ij^}), U\ and U% are respectively, 
the on site and nearest neighbor Coulomb interaction. In 
Fig. 2 we show a three dimensional ring (shaded region) 
with a flux 4> through the center of the ring that can 
cause polarization. At a particular point r on the ring 
we can bring an STM tip at a voltage V to cause further 
polarization at that position r (or i) while another STM 
tip can measure the local potential V(r) (or Vi) at r (or 
i). The polarization charge in a segment of the ring can 
be measured by a cylinder around the ring by looking at 
the induced charge on this cylinder (unshaded contour in 
Fig. 2). 



III. ANALYTICAL DERIVATION 

We will outline here all the mathematical steps re- 
quired to describe the polarization of a quantum system 
in terms of electrostatic quantum capacitance. When as- 
sumptions are used, we will give numerical verification 
and also cite appropriate earlier works. Suppose the po- 
tential V(r) at a point r is changed infinitesimally giving 



rise to a delta potential term in the Hamiltonian, of the 
form dV ext (r) = J2 n nS(r — r n ), where r„ is the coordi- 
nate of the nth electron. The Kohn-Hohenberg theorem 
[l|| states that the energy is an unique functional of the 
local potential. But we just use first order perturbation 
correction to the energy which is in terms of the applied 
potential only and no self consistency is required. So the 
increase in energy of the system can be expanded as 
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As k — i- 0, then 

dE 



dV(r] 



= Q{r) 



(3) 



where 

Q{r) = M J d 3 r 2 ...d 3 rM^*(r,r 2 ,...r M )^(r,r2,....,rM) 

(4) 

is the charge at r. We verify this numerically for all the 
geometries shown in Fig. 1. A plot is shown in Fig. 
3 for a disordered ring whose site energies ti vary from 
— 0.5t to 0.5t. Other parameters are explained in figure 
caption. We have used a single disorder configuration 
as the agreement is equally same for all other configura- 
tions. The figure shows the correctness of Eq. 3. We 
stress that we use exact diagonalization using Lanczos 
algorithm to determine E and Qi and hence this is a 
numerical verification of Eq. 3. 

Therefore, we can define a Lindhard function n(r). 
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The last step follows from Eq. 3. For the geometry in 
Fig. 1(a) or 1(b) or 1(d) or 1(e), we can change the 
magnetic flux (j) through the center. This will immedi- 
ately cause a redistribution of electronic charge in every 
site of the system. This is shown in Fig. 4 for the site 
numbered 8 in a 11 site ring. There is no qualitative 
difference between Fig. 1(d) and 1(e) even when the 
stub is weakly coupled to the ring. Also consider a case 
when an external voltage dVf xt is applied at site i (in 
the discrete model one can change ti infinitesimally as Vi 
and ti are linearly related). Obviously one can change 
both, (f> and V^ xt simultaneously. Charge at site i, Qi 
will change. Due to electron-electron interaction the lo- 
cal potential Vi, will change. Such a potential change will 
in turn have a feedback effect on charge displacement to 
determine dVi . This feedback effect is a purely quantum 
effect as this feedback occurs because the initial charge 
displacement (in this case induced by externally changed 
flux or potential) can change the quantum states or eigen- 
energies of the system. The externally applied flux and 
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potential also work indirectly by affecting the quantum 
states of the system. In other words external changes 
give rise to electron displacements for which potential at 
site i (or j) change, and this in turn gives rise to further 
charge displacements that is self consistently determined 
by Coulomb interactions and Fermi statistics. dV. (or 
dVj ) is this self consistently determined increment in po- 
tential at site i (or site j) and dQi is the self consistently 
determined charge at site i. We can begin by writing for 
the discrete system 



dQt = 



dQ t 



aV, 
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Although Qi is a functional of Vi, one can change the 
potential infinitesimally at a particular point without 
changing the potential at any other point to define a 
partial derivative and this is in fact done to arrive at 
the concept of functional derivative 19]. Not to men- 
tion, that in the standard definition for total derivatives, 
-j^- is change in Qi due to an infinitesimal test change in 
Vi "(i.e., Vf xt ), but dQ, on the LHS and dV and dVj on 
the RHS are actual changes which in these systems are 
determined self consistently. By the virtue of the fact 
that we are considering a sum over j makes our treat- 
ment valid for any geometry and any dimension, where 
the sum over j will run over all the sites making the sys- 
tem. For the continuous system Eq. 6 becomes, where 
the sum is replaced by an integration, 



dQ(r) = 



dQ(r) ,, , dQ(r) 



00 



dV(r) 



dQ{r) dV(r')d 3 r> 
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Here partial derivative with respect to V(r) means we are 
changing the potential in the region r to r + dr infinites- 
imally [19]. In the following we will argue analytically 
and numerically that 



E dQj 
dV~ 



Or by replacing the sum by integration, 
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While LHS is a sum of partial derivatives, the RHS is a 
functional derivative. Integration or sum of all (but one 
point) partial derivatives on the LHS is approximately 
a total derivative. Had this one point (r = r') been in- 
cluded it would have been an exact total derivative with 
respect to e nerg y (electronic charge times potential be- 
ing energy) [201 ] . So in the above approximation an en- 
ergy derivative is being replaced by a functional deriva- 
tive with respect to local potential. This approximation 
is known in other context like deriving the semi-classical 
limit of Friedel sum rule [2p| ] . That means this approx- 
imation is valid in the semi classical regime. Numerical 
verification of the approximate equality in Eq. 8 above 



is shown in Fig. 5. By doing exact numerical diagonal- 
ization it is difficult to go to a truly semi-classical limit 
as that will require us to take a large ring. However, the 
approximation in Eqs. 8 and 9 are known in context of 
deriving the semi-classical limit of Friedel sum rule [2(| ■ 
If this approximation holds then one can relate induced 
voltage and polarization charge through quantum capac- 
itance as shown below. Therefore, from Eq. 6 and 8, 



dQi 



dQ, 



. dQ, 
dVi 



SV, 



(10) 



Similarly, from Eq. 7 and 9, 



(ii) 

We have to assume that dVj is independent of j as further 
explained below. Now it follows from charge conservation 
that 



dQi \ - 6Qi_ 
'h 1 SVj 



dV 



(12) 



The RHS is the net change in Qi due to an infinitesimal 
functional increase (or decrease) in the potential at all 
points except at i. That is equivalent to not changing 
the potential anywhere but decreasing (or increasing) the 
potential at i infinitesimally. Due to charge conservation 
the change in charge at i will be the same in both cases. 

Coulomb repulsion tends to distribute charge uni- 
formly in a system. On the other hand quantum inter- 
ference effect tends to give rise to un-even distribution 
of charge. Assuming that in the semi-classical regime, 
Coulomb interactions dominate over quantum interfer- 
ence effects and distributes the charge uniformly, dVj be- 
comes independent of j and we denote it as dV res t . Nu- 
merical calculations for small size quantum systems show 
that for a wide range of parameter space the charge dis- 
tribution is uniform [2ll [22j. Only in very low density 
regime, quantum interference effects dominate and the 
charge density breaks up into crests and troughs [2l[ 1 2 21 ] . 
Any measurement process may not be able to resolve 
these crests and troughs and may show an average value 
for the local potential implying dVj can be taken to be 
independent of j. Transverse variation can be mapped 
to an effective variation in the longitudinal direction [23j . 
So one can write to a linear order, 



dQ, = C(dVi - dV rest ) 



(13) 



where C is the definition of geometric or classical capac- 
itance. We know that when we ignore quantum inter- 
ference effects (i.e., large systems without boundary and 
impurity effects averaged out or treated in random phase 
approximation) we can always get such a linear regime. 
Substituting Eq. 12 and Eq. 13 in Eq. 10 and simplifying 
one gets that 
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(r AdVj -dV rest ) = dQ l = 9Qj «9V^est , 

1 ^ <ty 90 ~ OVrest 8(f) [ ' 

Now since V res t can be changed by changing V. or <p, one 
can write V^ es t(</>, Vi). Therefore, 

„ r dV rest dV rest 
dVrest = q 7 »<? H ^77~dVi 

dip oVi 

Since, the region indexed z is very small compared to the 
rest of the system, a ^ 3t — ?-0 Therefore, 



dVrest dV res t 



(15) 



From Eq. (14) it follows that 

, \ (dVi — dVrest) dVrest dV re st 

iC + ??) T4> = = 

Multiplying both sides of the above equation by C we 
get, 



C 



(dVi - dVrest) C V dV r , 



C + 17 



or on using Eq. 13 



Crj 

dQi = 77~~ dVrest 

C +7] 



or 



dQ t . = —dQ r 



CeffdVr, 



where, C e// = That i 
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(16) 



(17) 



(18) 



When we define capacitance we do not consider the sign 
of the charge. Normally one plate of the capacitor has 
charge +Q and the other has charge — Q, wherein we 
write Q = CV . Hence if we want to include quantum 
effects then only in the semi-classical regime we find it 
possible to describe polarization in terms of capacitance. 

A quantum capacitor of capacitance r\ in series with 
the classical capacitance determines the effective capac- 
itance of the system. The characteristic potential (or 
the potential difference between two parts of the ring) is 
determined by this effective variable. The AC response 
of the ring is also determined by this effective variable 
along with the inductance. For a purely capacitive re- 
sponse, I(t) = dQi/dt or dl u = —iujCeffdV^ [l3j|. Such 
an effective variable will exist only if assumptions given 
in Eq. 3, Eq. 8, and Eq. 15 are valid. Although quantum 
capacitance was introduced first by Serge Luryi [24|, he 



introduced it on very general grounds and the above rela- 
tion was not obtained. The above relation was obtained 
for a closed system only in the frame work of single par- 
ticle two level system We have derived it generally 
for any arbitrary system including many body effects and 
outlined the assumptions so that such an effective vari- 
able will exist. 

In Figs. 6 and 7 we have made a comparison of quan- 
tum capacitance at site 8, i.e., rjs in the non- interacting 
system and the interacting system to show that r/ is 
a good parameter to effectively capture the effect of 
Coulomb repulsion. Here again we are considering the 
11 site disordered ring considered in Figs. 3 and 4. We 
have applied a small external potential at the 8th site to 
evaluate dQs/dVs = rjs- Fig. 6 is for the non-interacting 
case i.e., U\ = and U2 — although the electrons still 
interact through Fermi statistics. In Fig 7 we have made 
U\ = 2 with all other parameters remaining the same. 
There is a large qualitative as well as quantitative differ- 
ence between the two figures which shows the importance 
of including Coulomb interaction and many body effects 
in defining capacitance. dQi — —dQ re st can be measured 
as outlined in Fig. 2. = —77 can also be measured 
as outlined in Fig. 2. Ceometrical capacitance C is in- 
dependent of finite size or quantum interference effects 
and is known for a given sample from its bulk properties. 
So one can determine C e ff- Thus dVi and dV res t can be 
known for any applied external potentials as C e f / is the 
single parameter that determines this. 



IV. CONCLUSION 

In conclusion, in the semi-classical regime we prove po- 
larization charge and induced potential of a mesoscopic 
isolated sample are related by an effective capacitance 
C e f f ■ Effective capacitance can be decoupled as a linear 
combination of classical capacitance and quantum capac- 
itance. The quantum capacitance is given by the Lind- 
hard function. In this regime, we can design quantum 
circuits in terms of this parameter C e ff just as classical 
circuits are built in terms of parameters like resistance, 
capacitance and inductance. While in earlier works , Eq. 
18 was derived for single particle theory for a two level 
system, we have started from the principles of many body 
physics and derived Eq. 18 for any arbitrary geometry in 
any dimension. We have shown the approximations nec- 
essary to get this in the framework of fully interacting 
fermions. All these approximations are justified in the 
semi-classical limit. So Eq. 18 can provide a simple way 
to understand DC and AC response of quantum finite 
sized many body electronic systems in the semi-classical 
limit, in terms of an effective variable. 

The authors acknowledge useful discussions with Dr. 
P. Koskinen and Dr. S. Gupta. One of us (PSD) also 
thanks DST for funding this research. 
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Fig. 1. Schematic diagrams of some mesoscopic ge- 
ometries used in this work as examples and also for 
numerical verifications. Fig. 1(a) represents a one- 
dimensional ring pierced by a magnetic flux <f>. Fig. 1 (b) 
represents a one-dimensional ring to which a quantum 
wire or stub is attached. The ring is again pierced by a 
magnetic flux </>. Fig. 1(c) represents a two-dimensional 
square geometry Once again a magnetic field can be ap- 
plied perpendicular to the plane of the geometry. Figs. 
1(d), 1(e) and 1(f) are discrete versions of those in 1(a), 
1(b) and 1(c) respectively, that can be described by a gen- 
eralized Hubbard Hamiltonian and useful for numerical 
verifications. The dots represent sites. Nearest neighbor 
sites are marked % and j. 

Fig. 2. A three-dimensional mesoscopic ring pierced by 
a magnetic flux <f). The ring can be polarized by a voltage 
probe (without making a contact) at site r at a given 
voltage V. The flux can also polarize the ring. Another 



voltage probe (without making contact) whose voltage is 
allowed to vary can measure the voltage at the site r. 
A solenoid around the ring, as shown in the figure, can 
measure the induced charge in a segment of the ring due 
to polarization. The idea of quantum capacitance is valid 
only when the polarized charge is uniformly distributed 
in the rest of the ring, apart from the region at r. 

Fig. 3. The figure shows that for the geometries shown 
in Fig. 1, Eq. 3 is valid. In this figure we have used the 
geometry in Fig. 1(d). It consists of 11 sites (N = 11), 
with 4 spin up electrons and 4 spin down electrons. The 
on-site Hubbard U\ = 2, the nearest neighbor Hubbard 
1/2 = 1- The hopping parameter t = 1. The solid line 
is the charge on the 6th site as a function of the flux in 
units of Zq which is just electronic charge taken to be 1. 
The dashed line is j^j- in units of Zq = electronic charge 
taken as 1. is the ground state eigen-energy of the 
many body system found by exact diagonalization using 
Lanczos algorithm. The dimension of the matrix being 
of the order 10 5 AT10 5 . Here O = hc/e. We have also 
verified Eq. 3 for the geometries in Fig. 1(e) and 1(f). 

Fig. 4. The fig shows that a disordered mesoscopic 
ring can be polarized by an Aharonov-Bohm flux alone. 
We have used the geometry of Fig 1 (d) here. Parameters 
used are shown in the inset. We have used only one disor- 
der configuration for this figure with site energies varying 
from — M to M, but we have also checked for other dis- 
order configurations. The graphs for other geometries is 
qualitatively similar and so not shown here. Qg is the 
charge density at the 8th site, in units of Z which is 
electronic charge. At zero flux we expect the system to 
be neutral. As the flux changes, strong dispersion of Qg 
suggests polarization of the system wherein the positive 
charge in the system can be taken to be uniform and 
independent of flux. Here </>o = hc/e. 

Fig. 5. Here we are considering a ring in the semi clas- 
sical limit. That is the potential in the ring varies very 
slowly compared to de-Broglie wavelength. In particular 
we have taken a 11 site ring with a single defect, t\ = 0.3 
and rest of the site energies are 0. The solid line is 

and the dashed line is - fy 3 "- Both quantities are in 
units of electronic charge taken as 1. Here Ui = 2 and 
U2 = 1 with 4 spin up and 4 spin down electrons in the 
ring. Here 4>o — hc/e. 

Fig. 6. The figure shows a plot of rjg — dQg/dVg as 
a function of flux. Here rfl = e/t where e is electronic 
charge. <j>o = hc/e. This is for an eleven site disordered 
ring with site energies varying from — .ht to +.5i. In 
this case U\ = and U2 = 0. Which means electrons 
are interacting only through Fermi statistics. Here <po = 
hc/e. 

Fig. 7. The figure shows a plot of ijg = dQg/dVg as 
a function of flux. Here 77° = e/t where e is electronic 
charge. <f> — hc/e. This is for an eleven site disordered 
ring with site energies varying from — .5t to +.5i. In this 



case Ui — 2 and U2 = 0. Here <po = hc/e. 
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